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We present a measurement of the tt production cross section multiplied by the branching ratio to 
tau lepton decaying semihadronically (th) plus jets, cr(pp — > tt + X) ■ BR(tt — > ru +jets), at a center 
of mass energy ^fs = 1.96 TeV using 1 fb^ 1 of integrated luminosity collected with the DO detector. 
Assuming a top quark mass of 170 GeV, we measure a t j ■ BR^+j = 0.60lo'22 (stat) Vo'14 (syst) ± 
0.04 (lumi) pb. In addition, we extract the tt production cross section using the tt — > th + jets 
topology, with the result a t t = 6.9 t.1'2 (stat) +g'* (syst) ± 0.4 (lumi) pb. These findings are in good 
agreement with standard model predictions and measurements performed using other top quark 
decay channels. 

PACS numbers: 13.85.Lg, 13.85.Ni, 13.85.Qk, 14.65.Ha 



The decay t — ¥ Wb — ¥ rv T b provides a unique labora- 
tory to investigate the properties of the third generation 
fermions — the top (t) and bottom (b) quarks, the tau 
lepton (r), and the tau neutrino (v T ) — in a single pro- 
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cess. In the standard model (SM) , the t quark branching 
ratio (BR) to a W boson and a b quark is w 100%, and 
the final state is determined by the SM BR of the W 
boson. Since the t is the heaviest quark and the r the 
heaviest lepton, any non-SM mass- or flavor-dependent 
couplings could change the t quark decay rate into final 
states with r leptons. Therefore, any deviation in the 
BR of t — > Tv T b from that predicted by the SM can be 
an indication of non-SM physics. For example, in the 
Type 2 two-Higgs doublet model [l[ , such as required by 
the minimal supersymmetric standard model [2j, the t 
quark can have a significant BR to a charged Higgs bo- 
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son (-ff 1 * 1 ) and a b quark if m H ± < m t — nib- For large 
values of tan j3, the ratio of the vacuum expectation val- 
ues of the two-Higgs doublets, the charged Higgs boson 
preferentially decays to tu t , thereby increasing the BR 
of t — > Tv T b relative to the SM expectation and leading 
to a larger measured a(pp — » tt + X) ■ BR(ti — > t + jets) 
compared to the value expected from SM assumptions for 
the BRs and the production cross section [3J45J]. Other 
possible non-SM processes that can enhance the t quark 
to r lepton BR are .R-parity violating decays of the t 
quark in supersymmetric models [6J and new Z' bosons 
with nonuniversal couplings 

In this article, we present the first measurement of tt 
production in the r + jets final state using a data sam- 
ple corresponding to an integrated luminosity of 1 fb _1 
collected with the DO detector [§] at the Fermilab Teva- 
tron pp Collider operating at a center of mass energy 
y/s = 1.96 TeV. This measurement uses semihadronic r 
lepton decays, with BR ps 65%, as secondary electrons 
and muons from r lepton decays are difficult to distin- 
guish from primary electrons and muons resulting from 
W decays. Previous measurements of tt production using 
r leptons in the final state have been performed by the 
DO [sl and CDF [l(| collaborations in the Th+t channel, 
where t/j represents semihadronic r lepton decay modes 
and t represents either an electron or a muon. 

We apply the following preselection requirements: 
events must satisfy a multijet trigger requiring at least 
four jets; this is the same trigger used in the tt cross sec- 
tion measurement in the all-hadronic decay mode 
Reconstructed events are required to have missing trans- 
verse energy E T > 15 GeV and E T significance > 3, 
where the $t significance is a measure of the likeli- 
hood that the $t arises from physical sources rather 
than fluctuations in the measurement of the energies of 
the physics objects (jets, muons, electrons and unclus- 
tered energy) [l2j. Each event must also have at least 
four reconstructed jets with pseudorapidity |ry| < 2.5 
and transverse momentum pr > 15 GeV using an it- 



erative jet cone algorithm 
V(A 7? ) 2 + (A0) 2 = 0.5 [l4| 



with a cone size A1Z — 
The jet energies are cor- 



rected for the energy response of the calorimeter, the 
cone size, multiple pp interactions, event pile-up, and 



calorimeter noise [15J. At least one jet is required to 
have pr > 35 GeV, and at least two jets are required to 
have pr > 25 GeV. Each event is also required to have at 
least one Th candidate with pt > 10 GeV, \p\ < 2.5, and 
tau neural network output, NN T > 0.3 [l6j. Finally, to 
ensure this analysis is statistically independent of other 
DO tt cross section measurements so that it can be in- 
cluded in a combined cross section measurement, events 
satisfying the requirements of the tt — > e(fi) + jets chan- 
nel 17j, which include an isolated electron (muon) with 
Pt > 20 GeV, are rejected, as are events satisfying the 
requirements of the tt cross section measurement in the 
all-hadronic channel 11 1. 



A semihadronic r lepton candidate is a calorimeter 
cluster of cone size AJZ = 0.5 that includes any sub- 
clusters that might be present with E > 800 MeV con- 
structed from cells in the electromagnetic (EM) sec- 
tion of the calorimeter and the associated tracks with 
Pt > 1-5 GeV in a cone A1Z — 0.3 contained within the 
calorimeter cluster. These r candidates are classified ac- 
cording to one of three types based on the number of 
tracks and activity in the EM calorimeter, motivated by 
the semihadronic r lepton decays: (1) r — > n v T , (2) 
r ± -> ti^ttV, (3) r ± ->■ 7r ± 7r ± 7r = F( 7 r°) I / T . We define the 
three tau-types as follows: a single track with no EM 
subclusters (tau- type 1); a single track and > 1 EM sub- 
clusters (tau- type 2); and at least two tracks and > EM 
subclusters (tau- type 3). 

To further reduce the number of quark and gluon jets 
reconstructed as r leptons, we train separate neural net- 
works for each lepton decay type to improve the dis- 
crimination of t lepton candidates from the jet back- 
ground. The input variables to NN T are chosen to be 
minimally dependent on the t lepton energy and to ex- 
ploit the low track multiplicity and the narrow width of 
the calorimeter cluster produced by r leptons decaying 
scmihadronically, the low mass of the r lepton, and the 
differences in longitudinal and transverse shower shapes 
between r leptons and jets [l6| . A total of 12 NN T in- 
put variable are used to characterize the presence and 
properties of Th leptons, with seven of these variables in 
common for all three tau-types. The 12 variables are clas- 
sified as follows: isolation variables, shower shape vari- 
ables, and correlation variables between the calorimeter 
cluster and the associated charged particle tracks. Each 
NN T is trained on Z ^ t + t~ Monte Carlo (MC) events 
for signal and jets from data, where a jet and a noniso- 
lated muon are back-to-back in </>, for background. These 
are the same training samples used in Ref (l8l |. 

To measure the number of tt — > r^+jets signal events in 
data, the physics and instrumental backgrounds must be 
determined. The main physics backgrounds are H^ + jets 
events, where the W boson decays to a r lepton, and to a 
smaller extent Z + jets events, where the Z boson decays 
to a pair of r leptons with one misidentified as a jet and 
the 1£t is due to the neutrinos from the decays of the r 
leptons. The main instrumental background is multijet 
production where a jet is misidentified as a t lepton and 
the energy is mismeasured leading to a net $t- 

The preselection efficiencies and SM BRs for tt to fi- 
nal states with leptons [l{| are given in Table Q] These, 
as well as the final efficiencies, are calculated using a 
MC simulation of the experiment. The tt signal with 
leptons in the final state and W(Z) + jets background 
are simulated using the ALPGEN 1 . 2 [2(| matrix element 
generator assuming a t quark mass of 170 GeV and us- 
ing the CTEQ6L1 [2l[ parton distribution function set. 
These events are then processed through pythia 6.2 [22| 
to simulate parton showering, fragmentation, hadroniza- 
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TABLE I: A summary of the SM BRs of the various it sub- 
processes and the preselection efficiencies, where the uncer- 
tainties are derived from MC statistics. The leptonic r lepton 
decays are included in the e and fi channels, and represents 
an e, (i or t lepton. 







BR (%) 


^preselection (%) 


tt - 


+ Th + jets 


9.75 


40.5±0.2 


ti- 


■> e + jets 


17.7 


17.0±0.2 


tt- 


+ fi -+- jets 


17.6 


ll.liO.l 


it- 


■> Z + P + jets 


11.1 


4.04±0.03 



tion, and decays of short lived particles, except for b 
hadrons and r leptons. evtgen [23| is used to model 
the decays of b hadrons, while r leptons are decayed us- 
ing TAUOLA To avoid double counting final states 
generated by the leading-order parton-level calculation 
of ALPGEN and the parton-level shower evolution of 
PYTHIA, a matching algorithm is used [25[. The gen- 
erated events are then processed through the GEANT- 
based [26} simulation of the DO detector providing track- 
ing hits, calorimeter cell energies and muon hit informa- 
tion. The same reconstruction algorithm is applied to 
data and simulated events. 

The preselected data sample is used to extract the sig- 
nal and to study the multijet background after additional 
selection criteria are applied. To extract the signal sam- 
ple, we require NN T > 0.95. The selected events are 
then separated on the basis of tau-type according to the 
r lepton candidate with the highest value of NN T . This 
is done primarily to separate tau-type 3 events from the 
tau-type 1 and 2 events, since the former has a much 
higher misidentification rate and thus result in larger un- 
certainties on the tt cross section. In addition, we re- 
quire that each event have at least one identified b jet us- 
ing the 6-tag neural network {NNb) with the requirement 
NN b > 0.775. The NN b uses nine input variables that 
characterize the presence and properties of secondary ver- 
tices and track impact parameters within the jet [27j | . 
The efficiencies of these selections are shown in Table [TTJ 

The expected fraction of it events in the signal sam- 
ple is ~ 15% for tau-type 1 and 2, and 3% for tau- 
type 3 assuming a t i — 6.9 pb as measured in this analy- 
sis. In addition, the signal sample contains W(Z) + jets 
and multijet background events that must be subtracted. 
The W(Z) +jets contamination is determined using MC 
events, while the multijet background is determined from 
data. We start with the preselected sample and apply 
a loose t lepton veto, NN T < 0.9. Using MC events, 
we expect that the resulting sample contains < 2% 
tt — > Th + jets events and < 3% W{Z) + jets events, 
and therefore provides a good representation of the mul- 
tijet background. To further improve the modeling, the 
W(Z) + jets expectation is subtracted from the multijet 
background data sample. 



The numbers of signal and background events are ex- 
tracted from the final selected sample using a neural net- 
work (NN s b) event discriminant with the following input 
variables: (1) the scalar sum of the px of all jets and the 
t lepton candidate in the event; (2) the aplanarity [28j |; 
(3) the |?t significance; (4) the invariant mass of all jets 
and the r lepton candidate in the event; and (5) a x 2 
representing how well the 2 and 3 jet invariant masses 
agree with values expected for hadronic t quark decays, 
X 2 = (M 3 jet - m t )/cr 2 + (M 2 jet - mw) 2 /a w , with M 2 j et 
(M 3 j et ) being the 2 (3) jet invariant mass, m, t = 170 GeV, 
a t = 45 GeV and m w = 80 GeV, a w = 10 GeV are the 
mass and its resolution in the all-hadronic final state for 
the t quark and W boson, respectively. The jet combi- 
nation minimizing the x 2 is used. The NN s b is trained 
using a generated it — > Th + jets MC sample for signal 
and half the multijet data sample for background. 

We apply the trained NN s b to the signal data sample, 
the remaining half of the multijet sample, a tt MC sample 
with leptons in the final state that is independent of the 
NN s b training sample, and a W(Z) + jets MC sample. 
The application of NN s b on the multijet and MC samples 
is used to generate templates, as shown in Fig.Q] that are 
used to determine the fraction of it and multijet events 
using a negative log-likelihood fit. The normalization of 
the W(Z) + jets MC sample is derived by scaling the 
W(Z) transverse (dilepton) mass distribution to data. 
The normalization for it — ¥ e(/i) + jets is fixed to the 
theoretical cross section || and BRs. 

The number of it — >• Th + jets events extracted from 
the fit to data are 25.1±iJ§ (stat) and 18.01JJJ (stat) 
for channels with tau-types 1 and 2 together, and with 
tau-type 3, respectively. The fitted numbers of the 
multijet background events are 336.4±^J (stat) and 
1083.2tio 3 (stat), for the two channels, respectively. The 
numbers of it events are comparable to the expected val- 
ues given in Table HTT1 

To minimize the statistical uncertainty of the measure- 
ment of a(pp —> ti + X) ■ BK(ti —> Th + jets), which we 
denote as <j t i ■ BR r(i+J -, we fit the entire NN s b output dis- 
tribution rather than counting events above a given value. 
The value of cr t j ■ BR Tji+:/ and the fraction of multijet 
background in the sample are obtained from a negative 
log-likelihood fit to the NN s b distributions for tau-types 
1 and 2 and tau-type 3, independently: 

L(a th Nl , N° bs ) - - log (n ^T e " Wl J • 

where N { = a ti x £\ x BR tt - (i) x C + Nhk g ,i is the 
expected number of events in the i th bin of the JViV s j, 
histogram for a given cr t j, with integrated luminosity £, 
number of background events A^bkg.i, and the efficiency 
(BR) for the j th leptonic channel eti(j) (BRttm), and 
iV° bs is the observed number of events in the i th bin. 
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TABLE II: The efficiencies for the tight r lepton candidate (NN T > 0.95) and fc-tagging selections for tau-type 1 and 2, and 
tau-type 3 channels. The uncertainties are based on MC statistics. 



Tau-types 1 and 2 
Trigger (%) 



Tau-types 1 and 2 Tau-types 1 and 2 Tau-type 3 Tau-type 3 Tau-type 3 
NN T > 0.95 (%) fa-tag (%) Trig ger (%) NN T > 0.95 (%) fa-tag (%) 



tt — > Th + jets 
ti — > e + jets 
ti — > fi -+- jets 

-> ; + r + jets 



j+1.7 
-0.1 
3+1-5 
J -0.1 

65.9±S;? 
-+i.i 

-0.2 



74.8Z 
69.' 



50.5Z 



23.7±0.3 
33.1±0.4 
3.8±0.1 
43.7±0.4 



h2.8 
-2.7 

60.3+2:? 
1+2.8 
-2.7 



58. 7Z 



60.2Z 



?+1.8 
-0.1 

-1.4 
-0.1 
-1.1 
-0.2 



73.6] 



66.1] 
50. 7] 1 



19.4±0.2 
8.1±0.2 
7.7±0.2 

20.6±0.3 



+2.8 
2.7 
+2.8 
2.7 
+2.8 

61-4±H 



59.9 
58.9 
59.0 



10 2 
c 

iS 10 




Data (414 events) 
tt -> x h +jets 

Other tt leptonic decays 
| | W/Z+jets 
] Multijet 



-0.2 



0.2 0.4 0.6 0.8 



NN. 



TABLE III: Expected event yields in the two analysis channels 
assuming the measured ti production cross section of 6.9 pb. 
The uncertainties are derived from MC statistics. 



'sb 





Tau-type 1 and 2 


Tau-type 3 


tt — » Th + jets 


27.6±0.4 


22.1±0.3 


ti — > e + jets 


26.3±0.4 


5.9±0.2 


ti — > jj, + jets 


2.0±0.1 


3.7±0.1 


ti -> l + l~ + jets 


4.1±0.1 


2.0±0.1 


Total tt — > leptons 


61.3±0.6 


34.4±0.4 


W+jets 


13.5±0.3 


5.9±0.2 


Z+jets 


3.4±0.4 


1.9±0.1 



© 

c 10 2 
> 

LU 
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b) 




D0, 


1 fb -1 




+ 


• 


Data (1115 events) 












tt -> x h +jets 




i. 








Other tt leptonic decays 












W/Z+jets 










. I , 


Multijet 






-0.2 





0.2 0.4 0.6 0.8 


1 



NN. 



'sb 



FIG. 1: The output of NN sb for a) the tau-type 1 and 2 
channel, b) the tau-type 3 channel. The % 2 P er degree of 
freedom between data and templates is 0.6 for a) and 0.5 for 
b). 



The measured value of <j t i • BR 



IS 



0.60t°;^ (stat) (syst) ± 0.04 (lumi) pb, 



where we combine the tau-type 1 and 2 measurement 
with the tau-type 3 measurement. Using the theoreti- 
cal cross section a t t — 8.061q;y3 pb for m t — 170 GeV 
from Ref. [Ej], we measure BR Th+J = 0.074±g;g|f which is 
consistent with the SM value given in Table HI 

Table IIVI summarizes the systematic uncertainties on 
a t i ■ BR Tfe+:r These are calculated by varying the source 
by plus and minus one standard deviation, and propa- 
gating the uncertainty to the final a t i • BR rh+ j. The jet 
energy corrections account for the effect of the jet energy 
scale and resolution. Jet identification takes account of 
the difference in the jet finding efficiency in data and 



MC. The 6-tagging entry accounts for the systematic un- 
certainties on its efficiency. The r lepton identification 
uncertainty is derived by fluctuating the value of each 
input variable within its statistical uncertainty and ob- 
serving its effect on the NN T output. The trigger cate- 
gory accounts for the uncertainty in the multijet trigger 
turn-on and also takes into account the possibility that a 
multijet event with a r lepton can have a different trig- 
ger turn-on. Multijet modeling accounts for the uncer- 
tainty of the multijet sample to model the ti —> Th + jets 
background and its limited statistics. The category MC 
modeling accounts for the VK+jets modeling, the uncer- 
tainty in the scale factor both for light flavor jets and 
heavy flavor jets, and the parton distribution function 
uncertainty. The ti cross section systematic uncertainty 
represents the effect of the normalization of the non-tau 
lepton tt background, which is normalized to the theoret- 
ical value of the cross section. In addition to the sources 
listed in Table HV1 there is a ±6.1% uncertainty in the 
luminosity measurement [2{|. 

In addition, we present the combined measurement of 
the production cross section for ti using all measured ti 
channels with leptons in the final state listed in Table Hill 
that satisfy the selection criteria described above. We 
repeat the negative log-likelihood fit for the number of ti 
signal and multijet background events fixing the ti BRs 
to their SM values, but this time fit for all ti channels 
arriving at 60. 5± 11.8 (stat) events and 24. 0± 11.4 (stat) 
events for channels with tau-types 1 and 2 and with tau- 
type 3 characteristics, respectively. The fitted multijet 
backgrounds in this case are 336.7 ± 11.8 (stat) events 
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TABLE IV: Systematic uncertainties on a t f ■ BR T(i +j (in pb) as measured for the tt — > th + jets channel. 



Source 


i7i+jets (types 1 and 2) 


T,j+jets (type 3) 


Combined 


Jet energy corrections 


-0.078 


+0.081 


-0.047 


+0.047 


-0.068 


+0.069 


Jet identification 


-0.019 


+0.019 


-0.012 


+0.012 


-0.016 


+0.016 


b tagging 


-0.074 


+0.084 


-0.035 


+0.041 


-0.060 


+0.068 


Tau identification 


-0.035 


+0.035 


-0.020 


+0.021 


-0.029 


+0.029 


Trigger 


-0.002 


+0.053 


-0.000 


+0.027 


-0.002 


+0.043 


Multijet modeling 


-0.090 


+0.090 


-0.169 


+0.169 


-0.083 


+0.083 


MC modeling 


-0.028 


+0.028 


-0.012 


+0.013 


-0.023 


+0.022 


tt cross section 


-0.064 


+0.068 


-0.029 


+0.030 


-0.052 


+0.055 


Total systematic uncertainty 


-0.16 


+0.15 


-0.18 


+0.15 


-0.14 


+0.15 



and 1083.2 ± 11.4 (stat) events, for the two channels, 
respectively. The production cross section is calculated 
using the negative log-likelihood defined in Eq. [T]for tau- 
types 1 and 2 and tau-type 3 separately. The two cross 
sections are then combined to give 



ait 



6.9 ti:» 



(stat) to* (syst) ± 0.4 (lumi) pb. 



To estimate the dependence on m t , we reevaluate the 
efficiencies and templates using m t = 175 GeV and find 

a tl = 6.3 tH (stat) ±±°;£ (syst) ± 0.4 (lumi) pb. 

In summary, we have performed a measurement of a t i ■ 
BR Tfi+ j = 0.60i ; 26 an d) usm g the theoretical tt pro- 
duction cross section, extracted BR Th+ j — 0.074^q;q27, 
which agrees with the SM expectation. In addition, we 
have performed a measurement of the pp — > tt + X 
production cross section, a t i — 6.9 t{.i pb, using the 
tt — y T}i + jets topology. The measurement is in agree- 
ment with the SM 3 5] and previous experimental mea- 
surements using other tt channels [3] at the Tevatron. 
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